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]Q '. Abstract 

The transmissivity of a one-dimensional random system that is periodic on 
average is studied. It is shown that the transmission coefficient for frequen- 

o 
o 



cies corresponding to a gap in the band structure of the average periodic 
system increases with increasing disorder while the disorder is weak enough. 



This property is shown to be universal, independent of the type of fluctua- 
te | 

tions causing the randomness. In the case of strong disorder the transmission 
coefficient for frequencies in allowed bands is found to be a non monotonic 
function of the strength of the disorder. An explanation for the latter behavior 
is provided. 

PACS numbers: 42.25.Bs, 71.55.Jv, 73.20.Dx 
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It is well known today [1] that in the propagation of a classical wave through a one- 
dimensional periodic structure of finite length the amplitude of the transmitted wave de- 
creases exponentially with increasing length of the system when the frequency of the wave 
is in a gap in the photonic band structure of the infinite lattice of the same period. If this 
periodic structure is now randomly disordered in such a way that it remains periodic on 
average, it might seem as if a wave whose frequency lies in the gap would be even more 
strongly attenuated, due to the additional random scattering. In this paper we show, on 
the basis of three different models, that in fact the opposite occurs: the transmissivity of 
the disordered system is larger than that of the periodic structure, when the frequency of 
a wave propagating through it falls in a band gap of the latter. A possible reason for this 
counterintuitive result will be presented. 

We consider a scalar wave U(z) that satisfies the one-dimensional Helmholtz equation 

U{z) = 0, (1) 

where e(z) is a random, periodic on average, function inside the disordered region z e (0, L), 
and is equal to eo for z ^ (0, L). In the case when a plane wave of a frequency uj is incident 
on the segment (0, L) from the left the solution of Eq. ([l]) for z ^ (0, L) can be written in 
the form 

( e ik z + r (L\ e -ik z f OTZ< 

U(z) = (2) 
[ t(L)e lkoZ for z > L, 

1 /2 

where fco = Co 00 1 c -, which defines the reflection and transmission coefficients r(L) and t(L), 
respectively. 

The first model of a disordered structure that we consider is an alternating array of 2N 
dielectric slabs of dielectric constants e\ and e<i in the region (0,L). The width of the ith 
slab is given by = a(l + Aj), where the A.; are independent random variables that are 
uniformly distributed in the interval (—A, A), where clearly < A < 1. This structure is 
imbedded in a homogeneous medium, and a wave of frequency u is incident normally on 
it from the left. The complex transmission coefficient, t, has been calculated as a function 
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of the length L of the disordered structure by means of a transfer matrix approach [[J, for 
two frequencies of the incident wave in the vicinity of the lower edge of the lowest frequency 
band gap in the band structure of the average periodic lattice. For the parameters assumed 
(ei = 2, e 2 = 3.5) the frequencies of the lower and higher frequency edges of this band gap 
are defined by u_a/c = 0.873 and u + a/c = 1.038, respectively. In Fig. [IJ we plot the large L 
limit of — (a/L)(lnT) = a/l(u) as a function of the disorder parameter A for two different 
frequencies uj of the incident wave. Here T = \t\ is the transmissivity of the structure, the 
angle brackets denote an average over 100 realizations of the structure, and I (oj) is the 
Liapunov exponent. 

From the results presented in Fig. [I] we see that l~ l (uj) decreases monotonically with 
increasing A for the frequency of the electromagnetic wave within the band gap, which in 
turn means that the transmittance of the structure increases with increasing disorder. It is as 
if in the presence of the disorder channels for the propagation of waves in this frequency range 
that are closed in the absence of the disorder open up due to the partial filling of the density 
of photonic states in the gap region by the tails of this density of states from the higher and 
lower frequency bands bordering the gap. However, for the frequency in the allowed band, 
we obtain an at first sight surprising nonmonotonic dependence of the Liapunov exponent 
on the strength of the disorder A. The initial increase of l~ l (uj) with increasing A is readily 
understood as the suppression of the transmission in this frequency range due to the multiple 
scattering of the incident wave caused by the disorder. The decrease in the value of / _1 (co>) 
for strong disorder can be explained by a rather simple argument that will be presented 
below. The merging of the curves 1~ 1 {uj) for all values of uj for values of A larger than about 
0.6 is a reflection of the fact that for such large values of A the scattering structure has 
lost all traces of its underlying periodicity, as a consequence of which the band structure 
has disappeared, and waves of all frequencies see the same disordered medium that is now 
homogeneous on average. 

The second model of a disordered structure that we consider is a system of 2N — 1 
dielectric slabs, each of thickness a, in which the dielectric constant q of ith slab is given by 
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€2i-i = e + 5ei,€2i = 1, where the Set are independent random variables that are uniformly 
distributed in the interval [—A, A]. In contrast with the first model of a disordered structure 
considered, in which the disorder parameter A was limited by physical considerations to the 
values < A < 1, there is no restriction on the (positive) values that A can take in 
the present model. In Fig. g we plot the dependence of the Liapunov exponent / -1 (a;) as a 
function of the disorder parameter A when the frequency of the incident wave is in the lowest 
frequency gap of the photonic band gap structure of the average periodic system, and when 
it is in the allowed band below it. If the frequency is in the band gap, we see that when A is 
increased from zero, l^iuj) initially decreases, indicating that the transmittance T = |t| 2 is 
increasing with increasing disorder. However, when A reaches a value of approximately 4.5, 
begins to increase monotonically with a further increase of A. In contrast, when the 
frequency of the incident wave is in the allowed band l^iuo) increases monotonically with 
increasing A. 

The third model of a periodic system with disorder considered in the present paper is a 
continuous one: its dielectric constant is given by 

e(z) = Acosqz + 5e(z), (3) 

where 8e(z) is a zero- mean, Gaussian random process, with the correlation function 

(6e(z)5e(z')) = 5 2 e W [-\z-z'\ 2 /f COI }. (4) 

For a numerical calculation of the transmission coefficient we used the exact invariant imbed- 
ding equations 

d ^ = ^k e(L)[e-^ L + r(L)e^ L } 2 (5a) 

^= t ~k e(L)t(L)[l+r(L)e™° L ] , (5b) 

subject to the initial conditions r(0) = 0, t(0) = 1. The results are presented in Fig. |3|. 

When Figs. [TJ, 0, and |3| are compared, it is apparent that if the disorder is weak enough, 
the transmissivity for frequencies corresponding to gaps exhibits a rather unusual universal 



feature independent of the type of fluctuations: the transmission coefficient T(uo) ~ e~ L l l{ ^ 

increases when disorder appears (A 7^ 0), and continues to increase when the fluctuations 

grow. For strong disorder (A ~ 1 in the first model, A > 4.5 in the second one) the behavior 

of T as a function of A depends drastically on the type of disorder. Along with the expected 

decreasing of T in the second and third models one can see a quite surprising increase of 

the transmissivity of the system with "positional" discrete disorder (first model) when the 

strength of disorder A increases. 

To understand the differences of the dependence of /~ 1 (c<j) on the parameter A in the two 

discrete models studied in the present work, let us consider the disordered segment (0, L) 

as a set of random scatterers, each of which is a dielectric slab with a random width (first 

model) or a random dielectric constant q (second model) that possesses random reflection 

and transmission amplitudes R ni & and T in d, respectively. In addition to the randomness 

in the individual scattering characteristics, we have in the first model random distances 

5i =a(l + A 2 j) between consecutive scatterers. Then the transfer matrix M for the entire 

system can be written as the product M = M 1 F 1 M 2 F 2 ■ ■ ■ M N _iF N _iM N , where 

(oi &\ » fe lkS * \ 
M, = and F t = (6) 

are the transfer matrices for the ith scattering slab and for the homogeneous space between 
consecutive scattering slabs, respectively. The effect of the interfaces in this system are 
incorporated into the definition of the matrices I Mi >. The 11-element of the product MM + 
is related to the transmissivity of the structure T = \t\ 2 by 

(MM + ) u = l-l. (7) 

In averaging the matrix MM + we can use the independence of the matrices with different 
subscripts. This means that in the first model we can average over the distance between the 
(N — l)th and Nth. scatterers independently of all the other scatterers: 

(\a N \ 2 + \(3 N \ 2 ) (2a N (3 N e i2kS ^) 
[2a* N (3* N e- i2k5 "-i) (\a N \ 2 + \(3 N \ 2 " 



r N-l 



1 + R N 

T, 



N 



2 (a N (3 N ) (e^- 1 



(8) 

T~ / 

If the disorder is strong enough (2 A ■ ka ^> 1), we can neglect the off-diagonal terms on 
the right hand side of Eq. @. By repeating this procedure for each of the random transfer 
matrices, we obtain 



W) = (^f) I ) • (9) 

From a comparison of Eqs. (|7]) and (|9]) we find that 

1 + 



T I 2 



1 + -Rind 1 A 1 



(10) 



For a long system (N — > oo) Eq. (10) becomes approximately 



which means that the localization length /(u;) is of the order of 



-l(lnT) = ln/i±^) . (12) 



l(u) N' ' \1-R ind , 

In obtaining this result we have approximated (InT) by — ln(l/T), and have ignored ques- 
tions of the self-averaging or non-self-averaging of the latter quantity to obtain a simple 
qualitative result. For weak individual scattering we obtain 



— - 2 (R ind ) . (13) 



l(u>) 

Therefore in the case of strong disorder the transmission coefficient of a disordered random 
system is determined completely by the mean value of the reflection coefficient of a single 
scatterer. In the first model Ri nd is a periodic function of the width d of the single slab, 
and the averaging in Eq. (|T3|) means just the integration of R md over d in the interval 
(a(l — A),a(l + A)). It is clear that for A ~ 1 (strong "positional" disorder) the increase 
of the interval of integration of a periodic function causes the decrease of (R- m d)- 
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In the second model (random e^) the increase of A enhances the strength of a single 
scatterer and obviously leads to an increase of (-Rind)- 

We have carried out independent numerical calculations of (-Rind) (Fig- |])- The results 
are consistent with the reasoning above, and provide an explanation for the behavior of 
for strong disorder (large A) depicted in Fig. [TJ. 

In conclusion, we have shown on the basis of three different models of a one-dimensional 
random structure that is periodic on average, that the transmissivity of the structure for 
waves with frequencies corresponding to a gap in the band structure it possesses in the 
absence of the randomness is increased by the randomness. In addition, we have given 
a qualitative explanation of the observed nonmonotonic dependence of the length on the 
disorder parameter A for waves whose frequencies are in an allowed band. 
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FIGURES 

FIG. 1. Liapunov exponent as a function of the disorder parameter A for the first model 
with "positional" disorder. Curve a corresponds to the wave frequency tua/c = 0.90 above the 
lower frequency edge of the first band gap (uj-a/c = 0.873), and b corresponds to the frequency 
uoa/c = 0.85 below the edge. 

FIG. 2. Liapunov exponent as a function of the disorder parameter A for the second model 
with random ej. Curve a corresponds to the wave frequency uja/c = 0.6 above the lower frequency 
edge of the first band gap, and b corresponds to the frequency uoa/c = 0.5 below the edge. For the 
parameters of the second model assumed (e = 9.0) the position of the lower frequency edge of the 
first band gap is at uo^a/c = 0.568. 

FIG. 3. Liapunov exponent as a function of the root mean square of the random function 
Se(z) in the third, continuous, model. The curves a and b correspond to wave frequencies above 
and below the lower frequency edge of the first band gap, respectively. 

FIG. 4. Plot of (.Rind}- Curves a and b were calculated for the frequencies that correspond to 
the frequencies of the curves a and b in Fig. 1. 
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